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fe ee ROP UerTON 


tiemecnmmOey Hadley and Whitin [Ref. 1] is a standard 
reference on mathematical inventory models. A particular 
class of models presented therein which has received wide 
acceptance is known as the lot-size, reorder point, or 
(O74), model, where r is defined as that inventory level at 
which a reorder for a quantity Q is initiated. Both Q and 
r are assumed to be continuous variables. The authors 
begin their treatment of stochastic customer demand by 
deriving two models referred to as the backorders case and 
the lost sales case, respectively. In the backorders case, 
all demand occurring when the system is out of stock is 
eventually filled, whereas in the lost sales case, this 
demand is lost forever. 

The models derived by Hadley and Whitin consist of 
long-run average annual variable cost functions, both 
labeled K(Q,r). Each function consists of a sum of three 
Ost expressions, the first representing a set-up or ordering 
cost, the second an inventory carrying cost, and the third 
a backorders or lost sales cost as appropriate. The authors 
assert that these functions are convex in the variables Q 
and r and their partial derivatives are equated to zero to 
determine unique and global minimum costs in each case. {In 
particular, they claim these functions are strictly convex 


in the case where lead time demand (that demand occurring 





between the reorder point and arrival of that order) is 
normally distributed. These assertions are unproven 
however, and are left to the reader as a series of exercises. 

The first exercise asks the reader to show that a 
Poseicwrady ceormewhich occurs in both models, say J(Q,r), 
is convex. It 1S obvious that the other terms are convex. 
Pemvemwrd then Follow that K(Q,r) is convex in Q and r from 
additivity. However, it was pointed out in 1964 by Veinott 
[Ref. 2] that J(Q,r) is not always convex. In 1969, Brooks and Lu 
(Ref. 3] considered the convexity in the backorders model 
and showed, for normally distributed lead time demand, that 
ie aero cOmyex fom alliyr > t, the mean lead time demand. 
Peso7 Sa. Ehe nonconvexity of the backorders cost function 
and the nonuniqueness of its solution was demonstrated by 
Minh [Ref. 4] for the case where lead time demand is normally 
eee loueca. sim Parclcular, Minh showed through counter- 
examples that K(Q,r) is not convex and, if a solution exists, 
its partial derivatives equated to zero have exactly two 
solutions vice one as claimed by Hadley and Whitin. Minh 
also succeeded in classifying one of the solutions as a 
relative minimum and the other as a saddle point. Further, 
Minh showed that feasible minimum solutions can be obtained 
where r < Hu. 

As Minh has suggested, his results have created doubts 
S-eemeocrey ana Waitin's claim that the lost sales cost 
Sieeronwis strictly convex for the case of normally dis- 


tributed lead time demand. Indeed, the only mathematical 





difference between the two models consists of an additional 
term in the inventory carrying cost expression of the lost 
Slesmeecoe ew tuleriom.. This term accounts for the additional 
Sectewineurerea DY Carrying aS inventory that portion of an 
arriving lot which would have been immediately dispatched 
to satisfy accumulated backorder demand, if such were per- 
freced. This 1s consistent with the authors' assumption 
that the expected number of lost sales is equivalent to the 
expected number of backorders, all other factors being 
equal. 

If, as will be shown in this thesis, K is not convex 
as claimed, then of course a different set of arguments must 
be provided for its minimization. In particular, the con- 
straint imposed by assuming a nonnegative reorder point 
will also have to be examined. Hadley and Whitin ignore 
this constraint in their development, leaving the impression 
that the solution is always to be found in the region 
r>0O, Q0> 0. Indeed, they are not clear as to their 
intent regarding r at all. In the beginning of their 
derivation, they require r to be positive; in the exer- 
cises dealing with convexity however, r is allowed to be 
nonnegative. Since r is assumed to be a continuous variable, 
there seems to be no rational basis for excluding r = 0 and 
it will be assumed herein that r is to be treated as a 
nonnegative variable. 

It is the purpose of this thesis to clarify the issue of 


convexity for the lost sales case and completely characterize 





its solution, including the effect of the constraint on 
r, assuming a normally distributed lead time demand. 
Section II will show that K is not convex and that the 
nonnegativity constraint on r precludes Hadley and Whitin's 
claim of existence in that the simultaneous solution to 
the partial derivations of K may be outside its domain of 
definition. Section III will present a series of lemmas 
and theorems in which the minimum solution is completely 
characterized. Section IV will briefly discuss these 
results. 

For a normally distributed lead time demand, Hadley and 
Whitin define the lost sales cost function as follows: 

Q 


KQx) = B+ rcS+r- yu] + (1e+ mA) (Ger) oY +0041 


The parameters are defined as follows: 


A is the set-up, or ordering cost. 

Ceplceene Unle Cost. 

MPEP tomnventory Carrying cost per unit per year. 
7 is the cost of a lost sale. 

A is the mean annual demand. 


G is the standard deviation of the normally distributed 
lead time random variable. 


Pio nenemnean O: that random variable. 


oz) is the standard normal density function. 
#z) is its complimentary cumulative distribution function, 
dae: D2) ee nd) (52) -Cii< . 





For ease in subsequent algebraic manipulation, K(Q,r) 


will be rewritten as 


EiGa, 2) = = + Ic(S + eZ + 77) ) + SO 
where Zz = — aoemez — Cidt{z) = zd(z)J]. In this notation, 


the set-up cost, inventory carrying cost, and lost sales 
costs are clearly identifiable. The domain of definition 
Botceigneds a EUneclon Of z and ©, is defined by the restric- 


Eons Z 7 = =, Oe >) te 


10 





PPC OUN TE REXAMPLES TO CONVEXITY 
AND EXISTENCE OF SOLUTION 
A function of two or more variables is strictly convex 
if and only if its Hessian matrix is everywhere positive 
definite. The Hessian matrix 1S positive definite if and 
only if its proper diagonal terms and determinant are 
Strictly positive. The Hessian for K(z,Q) is derived below. 


The Hessian matrix is defined as 


3°K 34x 
52° dZ90Q 
C= 
3°K °K 
j0Qez 30° 


Since K(z,Q) is continuously differentiable over the 








BeGlom 2 > - = oO= 0, 
sR ea 
aZ90Q 6002 


Poumot@eseduenc Gerivations, K will be written as 


Cz, 0). = =I + tn (2z)] + Ic (S Eee gete mz) | . 


Then, 


aK = TAy'(z) + IClo + 9° (2)1 


aa 








where 








fez y= -OO(Z) . 
Maus , 
dk _ _ TAod(z) a 
ae ae Geer oA) ey 

Hence, 

OK _ ; si, {OSS 8e 

a 0 implies DNGZ.= A= O*Te Fmt (iy 

°K Tea oz) 

OZ 

= OO(2Z) fore + oT A} 
Q 
°K | Mmicekz) a 3° 
dO0Z 2 dZ9Q 
O 
chee A 18€ 
wer i aoe + mu (z)J] + 7 
Q 

Then, 

aK _ . . e ay /ZALA + m(2*)d 

50 7 0 implies Q* = TC (2) 

ve 

oR = SSfa + m(z)] 

dQ Q 


The Hessian matrix 1s therefore: 


eZ 





Goce) (O1C + TA. TAG ® (z) 


0 9° 
2 = 
TAG O (z) 2h (ae (Zz) | 
oF Q° 


It is shown in the Appendix that the function n(z) 
MemseeiGtly positive. Thus the proper diagonal terms of 
Q are strictly positive and 2 is positive definite if and 
SMe vett det % {determinant of ©) is positive. The 
determinant of 2 1S given by 


det & = 2g) 20(2 Amen OLC tamu ee ra? (2) 
Q 


Define the expression contained within the outer brackets 


as) V¥(Z2,Q). Since ag Per Get. oP O alt and enly 2t 
Q 
2 2 
ene eee) teeter (zZ))] (OIC + mA] - ToAcd (z) > 0. 


Wmetethis result, a counterexample to the claim that K(z,Q) 
is strictly convex is available. The parameters utilized 
are from an example given by Hadley and Whitin where the 
lead time demand is normally distributed and are specified 
caw ee o00ene = 50, © = 0.20, A = 1600, n = 2000, o = 50, 
fae einen tos the policy r = 700 (2 = -1.0) and Q = 10, 


pez Oz 420), (2) = 0.8413 and so, 


1s 








Vez, ©) = 2(.2420) POU Or 200050) (22420 + 53413) | 


<[(20) (.20) (50) + (2000) (2600) | 


~ (2000)* (1600) (50) (.8413) 2 


= 1.7398214 x 10°) - 2.2649142 x 10t! < o. 


imc 2,0) 1S mOr Strictly convex throughout the region 
mae, © > 0 as claimed. 

As for a counterexample to the existence of the solu- 
tion, choose the following parameters: A = 25,000, 
ooo Omeie= 02.99, 71 = 10, A = 5, go = 10, and u = 30. 
For these parameters, the iteration procedure described 
Moaecemeey sand Whitin produces a solution (z2*,0*) = (-3.10,5.06). 


It may be verified by results to be subsequently derived 


that V(z*,Q*) is positive. It will be observed, however, 
Pete tne reorder point r corresponding to z* = -3.10 is 
me— eZ Oeteie— —1|,.0 which is outside the domain of definition 


for K(Q,r) as defined by Hadley and Whitin. Basically, 
the counterexample is a consequence of the fact that the 
solution (z*,0O*) is independent of u, the mean lead time 


demand. 
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Mite oleae ten ZATION OF THE MINIMUM SOLUTION: PART A 


It is convenient to separate out as lemmas certain 
isolated results that support the main results presented 
in this thesis. The first four of these lemmas are stated 
in this section and derived in the Appendix. They lead to 
an intermediate result stated as THEOREM 1, a theorem that 
demonstrates the uniqueness of the solution (z*,Q*) to 
equations (1) and (2) and establishes that solution to be 
an unconstrained relative minimum for kK. 

LeMmMAs 5 thzyu 7, also stated in this section and derived 
in the Appendix, lead to THEOREM 2, a theorem establishing 
PiIierorediieoolmts (7,0) where z ~*~ z*, K(z*,Q*) < K(z,Q). 
The main result will then be presented as THEOREM 3 in which 
the minimum solution is completely characterized. 

Poem ceesoreoarameters A, CC, I, wt, A, G, and up, all 
greater than zero, be given. Define the following functions 


mon —-~H < Z < @: 


Dew? 

Baez, ) = p AS (2) ~ Minyez oO 
Zic{l = o(z) ] 

G(z) = 3¢(z) + z - z0(z) 

H(z) = —?l2) - 
eo CZ) 


ibs 





LEMMA 1: 


LEMMA 2: 


LEMMA 3: 


LEMMA <¢:; 





(Z*,Q*) be a solution to equations 
andw 2) TANen, 
A= A(z*) 

> 
WetZ* Oo * } = 0 ie anc eo Ly a 

< 

> 

< 
Picomilimetlon G(Z)) LS Strictly positive. 
tice sunmectieonmehi(z) 1S strictly monotone 
decreasing. 

ree 

The equation H(z) = = has one and only 


one solution. 


HicmeimeetonoA(Z) has a unique Minimum at 


Zo where ZO is the unique solution to the 


(b) 





, Lee 
equation H(z) = — 

mana (2) =" O 

Zra 

Lim A(z) = © 

27-0 

A(Z.) < 0 
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PecOREM 1: 


mROUr : 


(d) There exists one and only one Zz, < Zo 
Such that A(z,) = 0. Moreover, if 


Za 7, AZ = 20% 


"ed 
There exists one and only one (z*,Q*) which is 

a Simultaneous solution to equations (1) and (2) 
and that solution determines a relative minimum 


tor, K. 


From LEMMA 1, let z* be chosen so that 
At(z*)}) =A > 0. Then z* < Z4 from LEMMA 4(d). 
Maeweeeme strict monotonicity of A(z) for 


EE Sas re Pewerenetollows Ghat’ z* 1s “unique. 


1’ 
fo waecseaen be selected according to either of 
eqiaerens (1) or (2). Then, from LEMMA 1, 
(z*,Q*) is a unique solution to these equations. 
To show that K(z*,0O*) is a relative minimum, 
NoOeeoenaz Z* < Z4 < Zo: Then, by LEMMA 3, 

H(z*) > 252. But LEMMA 1(b) establishes this 

aS a necessary and sufficient condition for 
positive definiteness of the Hessian of K at 
(z*,Q*). This positive definiteness establishes 


Eaoeew( 2 ©. jetssa crelative minimum of K and 


uniqueness has already been established. 


Sipe e Ie 


ey 





Iii. CHARACTERIZATION OF THE MINIMUM SOLUTION; PART B 


The results of the previous lemmas are depicted 
Graphically in Figure 1. The functions A(z) and H(z) and 
the lines A and aEe are shown for perspective in that they 
represent no particular set of parameters. It should be 
noted that z*, here shown negative, may be positive or 


negative. 


H(z) 


EUCZ ) 


Figure l 
It has been established in PART A that K(z*,Q*) 


constitutes a relative minimum and that (z*,Q*) is unique. 


iMmetsematt Le will be demonstrated that K(z*,Q*) does 


ies: 





constitute an unconstrained global minimum. However, 
fPameas Shown ae Section If that the constraint r > 0 
= =) may be binding. Therefore that constraint must 
be considered in the minimization of K. 

For the succeeding lemmas and theorems, it will be 


Semvenient to define the following functions for Q > 0, 


-~m~ < Z < @®, 


e TAD(Z) 
etme) =" O(Z) } 


zk 2A(A + tn (z) ] 
Qo (2) =\ Ic 


LEMMA 5: Let z be fixed. Then, K(z,Q,(z)) < K(z,Q) 


mer allv®@. 


pew os Lt 2 < Z*, Q, (Z) < Q, (2), anew Zee ez, 
Q(z) > Q, (2) 3 


HEitiA 7 seeet @ > 0 be fixed. Then,~1f z< Z1(Q), 


aK aK 
Ss Semana, 2f Zo 2 Zz, (Q), eras Oz 


Mis@mieves-. te 2 7 z*, then K(z*,Q*)< K{(z,Q) for 


siti Zzesanag sO > 0% 


ae, 





BROOE : 


Poueaeoe wht xed and vy = (cos 6, sin 6) be a 
Wms vector from (2,Q.,(z)) Mette. Girection Of 
(2 -7O-). ine directional derivative of % at 


(Z,Q.,(Z)) may be written as 


_ OK aK 
Pe lee co? 9 4G Sn 
= ok cos 86 
OZ 
; aK 
Since 30 = Q at (Z,Q.(2)). The proof now 


proceeds by cases. 


Gasoeyl: 2 < z*. Since Q. (2) iS Sstricely monovone 
decreasing, Q. (Z) > Q,(z*), and thus an < ss Air 
Hence, cos 86 > 0. Let Q1 = Q, (2) and Q. = Q, (Zz). 


By LEMMA 8, Q, < Q), in which case Z,(Q5) > 2, (Q)) 


Since Z4 1s strictly monotone decreasing. 


Since Z, (Q)) = Z, Z (Q,) 2 Zewnen. tmp les 

cS SSOeoymhe MAS. Hence, DV <Q) siehig 

Ze 

Pease. 2s. 2 > Ze. Since Q. (2) is strictly monotone 


decreasing, Q. (2) < Q,(2*), and thus 5 Mls) <5 ge 


Hence, cos § < 0. Let Q1 = Q, (2) and Q, = Q, (Zz). 
By LEMMA 8, Q, > Q,: in which case Zz, (Q.,) < Zz, (Q,) 
Since Z4 1s strictly monotone decreasing. Since 

, . dK 
Zz, (Q,) = Z, Z, (Q.) < zg which implies ap 2.0 LPYy 


LEMMA 9. Hence, DV < 10) St chau eer Ae 


20) 





MaHeOREM 3: 


PROOF: 


Pimeleecases, Dov < 0 sO that K is strictly 


K 
decreasing from any point (Z,Q.(Z)) toward 
OZ es) ee aus: K(z*,O*) < K(Z,Q.(Z)) and, by 


LEMMA 5, K(Z,Q. (Z)) ez 0) cen amy, 0 - = 0. 
CED: 


8 4s il = then K(z*,0*) is the global minimum 


Q 


Vameme Of KK. If z* < = =, then K(- =,Q,(~u/o) ) 


1s the global minimum value of K. 


SA = then the corresponding value of r 
1s nonnegative and r is in the domain of definition 
Ope eee AO ToPOnre 2, K(z*,O*) < K(z,O) for all 


@eezeeanad © > 0. Thus, K(z*,Q*) is a global 


Pine eo Oo and Z* > = =. 
Mia Zz =< = then the corresponding value of 


r is negative and is thus outside the domain of 
Goieeieemeneof K. Hence, K(z*,0O*) cannot be optimal. 
However, LEMMA 5 has established that, for any 
£1<ed Zz, K(z,Q.,(Z)) Sez, Cand, by THEOREM 2, 
K is strictly decreasing from any point (Z,Q.,(Z)) 


toward (z*,Q*). Hence, K({- =, QO, (- =) ) < K(z,Q,(Z)) 


mOtuealyy Zi). 7 = = Geese tne mrnimem value 
Spiieeoceurs On the boundary z= = = ang LOG 
z* < - 4 K is minimized at K(- a O.(- ine 
a’ Gace Oo 
Cao Oe 


Zi: 








iV-.) DISCUSSION 


The results of this thesis have resolved a long-standing 
Sesto Concerning the solution to the lost sales case of 
the Hadley-whitin (Q,r) model. Since the model is widely 
used by practitioners, it is an issue that needs resolution, 
Eermene seit ire argument given by the authors is based on the 
supposed convexity of the objective function. That not being 
the case, it is necessary to supply a different argument for 
the solution. Convexity aside, the boundary constraint can 
be binding and appears to have been completely ignored by 
the authors in their development. This thesis has clarified 
both these issues. 

Pevemameme mere Clarification of the results, this 
thesis exemplifies the fact that attention to the underlying 
mathematical detail and rigor cannot be ignored in modeling 
peewee o wen.) FOrtunately, in this case, the solution remains 
Poe yzevyalia but EOr a different set of reasons. This 
could only have been discovered by a close and detailed 
examination of the mathematics of the model. In the process 
of that examination, the role of the boundary condition has 
memeciicdenesel® and allows for complete characterization of 
tae solution. 

As for recommendations for further study, it should be 
pointed out that only the case of normally distributed lead 


time demand has been considered. There are, however, many 


eZ 





other assumptions that might be made concerning this demand. 
Among those of special interest to Navy applications: are 
those of the Poisson and negative binomial distributions. 

An investigation similar to that presented here for at least 
these two cases would seem to be called for and is a 


recommendation of this thesis. 


ase 





APPENDIX A 


PROOFS OF LEMMAS 


This Appendix contains the proofs of LEMMAs 1 thru 7 
inclusive. Prior to their proofs, it will be convenient 
to prove four auxiliary lemmas which are appealed to in 
the main argument. For clarity, the auxiliary lemmas will 


be numbered with the letter "A" following the number. 


LEMMA 1A: The function n(z) is strictly monotone decreasing 
for -™ < z < » and everywhere strictly positive 


over that region. 


PROOF: Meez y= Oil (Z) 
where 
Meeziee= O12) — 2ZO(Z). 
Now, 
Wane =e= 2012) =z (-o (Zz) ) = 942) 


=— Orez > <0 for alloz, 


ile Z ) OZ) ae. 0 POm a bie ss 


24 





Clearly, 





Pee Zien (ore) —ZO(7)) = 0 + og=—e +o, 
Zr Zr=— a 
And 
aenlZ eine OZ) —Z0'z)) = —Lim zo (Zz) 
Vda) Z>a Vda) 
= - Lim “s 
a On) 


Seo Mehdi so iain te tcesules 2n the indeterminant 


form =, L'Hopital's Rule is applied repeatedly. 

















: Z Bee eee Ll Eee ee Cz) 
mee el SS Tz) SS ey 
Zz Sz) Z 5 Zz 

6° (z) 

ee ee 2 Oz) 
= Lim —=Zzolz) 

7-0 
= - Lim aha OF 

Z 
7 
i Micmipizie— sOUlg) 1S a Strictiy convex decreasing 


PUMeELOnwamoroecning 0 as z++> and hence 


Wea s O £or all z: 


25 





LEMMA 2A: 


PROOF: 


LEMMA 3A: 


PROOF : 


LEMMA 4A: 


7 TA® (Zz) 
Q, (2) I@Ti ~ d(ay] iS strictly monotone 


decreasing and everywhere strictly positive. 


" 
' = IN se) OR) ->(z) 
Q CZ) ps _ > tO SS 
a IC} ry - 8(2))2 (= orlz)) 7 i 
Sa < 0 for all z. 


Ic{l - 6(2z)]° 


Q, (2) is everywhere strictly positive by inspection. 


Ore. 


Q. (2) = peste atai) is strictly monotone 


decreasing and everywhere strictly positive. 


Squaring Q. (2) andwcitterentiating lmplicitly, 


ATa (Zz) 


o, (ie Ge mae, az 


Q5 (2) = - 


Q. (2) is everywhere strictly positive by inspection. 


On aL. 


de Oe 
= oO A 
Z, (Q) f Orc + TA Weeser ice ly MoOnorone 


decreasing. 
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PROOF: Q, (2) 


Or 


-l OEE 
QIC + mA 


Since Zz (Q) is the functional inverse of Q, (2) 
and Q, (2) is strictly monotone decreasing, 


Zz, (Q) PS Sericely monotone decreasing. 


One). 


PEMMA 1: fa) SAe = Alzs). 


Proort (a). Equating the two expressions for Q* given by 


equations (1) and (2) yields 
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